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1 Abstract 

The concept of z-scaling reflecting the general features of internal particle sub- 
^r) . structure, constituent interaction and mechanism of particle formation at high-py 

is reviewed. Experimental data on inclusive cross sections obtained at U70, ISR, 
SpS, RHIC and Tevatron are used in the analysis. The properties of data z- 
presentation such as the energy independence, power law and ^-dependence are 
discussed. The properties of z-scaling are argued to be connected with the fun- 
damental symmetries such as self-similarity, locality and fractality. The use of 
z-scaling to search for new physics phenomena in hadron-hadron and hadron- 
. £h ', nucleus collisions is suggested. RHIC data used in our new analysis confirm 

z-scaling in pp collisions. High-pr particle spectra at LHC energies are pre- 
dicted. Violation of z-scaling characterized by the change of anomalous fractal 
dimension is considered as a new and complementary signature of new physics 
phenomena. 
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1 Introduction 



More prominent properties of particle are observed at high energy and transverse 
momentum [l]-[4]. The kinematic regime is used to perform calculations of physical 
quantities in the framework of perturbative Quantum Chromodynamics (QCD). Devi- 
ations of theoretical results obtained in a high-pr region using available experimental 
data are often considered as manifestations of new physics phenomena. At the same 
time we should note that nonperturbative effects are not well controlled by the the- 
ory. Therefore search for new scaling features of particle interaction in the region is of 
interest for development of the theory. 

The fundamental problem of high energy physics is the origin of mass, spin and 
charge of particles. The study of particle interactions over a wide kinematic range and 
especially at small scales is necessary to resolve the problem and understand underlying 
physics phenomena. One assumes that all types of interactions are unified at small 
scales. New ideas and theories such as extra dimensions, anisotropy and fractality of 
space-time, quark compositeness, theories of Grand Unification, Super Symmetry and 
Super Gravity are intensively developed. 

New feature (z-scaling) of high-pr particle production in hadron-hadron and hadron- 
nucleus collisions at high energies was established in [5]- [11]. The scaling function ip 
and scaling variable z are expressed via experimental quantities such as the inclu- 
sive cross section Ed 3 a/dp 3 and the multiplicity density of charged particles dN/drj. 
Data ^-presentation is found to reveal symmetry properties (energy independence, A- 
dependence, power law). The properties of tp at high z are assumed to be relevant to 
the structure of space-time at small scales [12, 13]. The function ip{z) is interpreted as 
the probability density to produce a particle with a formation length z. The concept 
of ^-scaling and the method of data analysis are developed for description of differ- 
ent particles (charged [6, 7, 14] and neutral [10, 11] hadrons, direct photons [8, 15], 
jets [9]) produced in high energy hadron-hadron and hadron-nucleus interactions. The 
proposed method is complementary to a method of direct calculations developed in 
the framework of QCD [16] and methods based on Monte Carlo generators [17]- [24]. 
Therefore we consider that the use of the method allow to reduce some theoretical 
uncertainties which are ambiguously estimated by theory. 

In the report the general concept of ^-scaling, the properties of data ^-presentation 
and some new results of data analysis are presented. The fundamental principles such 
as self-similarity, locality, fractality and scale-relativity are formulated and discussed 
in the framework of the z-scaling concept. Verification ^-scaling validity at RHIC and 
LHC is suggested. Violation ^-scaling is considered to be indication of new physics 
phenomena. 

2 Z-scaling 

In the section we discuss underlying ideas of ^-scaling. A general scheme of data 
^-presentation is described. Physical meaning of the introduced quantities is explained. 

2.1 Locality 

The idea of z-scaling is based on the assumptions [25] that gross feature of inclusive 
particle distribution of the process (1) at high energies can be described in terms of 



the corresponding kinematic characteristics 

M 1 + M 2 ^m l +X (1) 

of the constituent subprocess written in the symbolic form (2) 

(xiMi) + (x 2 M 2 ) -> mi + (xiMi + x 2 M 2 + m 2 ) (2) 

satisfying the condition 

(xiPi + x 2 P 2 - p) 2 = (xiMi + x 2 M 2 + m 2 ) 2 . (3) 

The equation is the expression of locality of hadron interaction at constituent level. 
The Xi and x 2 are fractions of the incoming momenta Pi and P 2 of the colliding objects 
with the masses Mi and M 2 . They determine the minimum energy, which is necessary 
for production of the secondary particle with the mass mi and the four-momentum p. 
The parameter m 2 is introduced to satisfy the internal conservation laws (for baryon 
number, isospin, strangeness, and so on). 

The equation (3) reflects minimum recoil mass hypothesis in the elementary sub- 
process. To connect kinematic and structural characteristics of the interaction, the 
quantity fl is introduced. It is chosen in the form 

Q(x 1 ,x 2 )=m(l-x 1 )^(l-x 2 ) 5 \ (4) 

where m is a mass constant and Si and S 2 are factors relating to the anomalous fractal 
dimensions of the colliding objects. The fractions x\ and x 2 are determined to maximize 
the value of Q(xi,x 2 ), simultaneously fulfilling the condition (3) 

dQ(xi,x 2 )/dxi\ X2=X2 ( Xl ) = 0. (5) 

The fractions x\ and x 2 are equal to unity along the phase space limit and cover the 
full phase space accessible at any energy. 



2.2 Self-similarity 

Self-similarity is a scale-invariant property connected with dropping of certain di- 
mensional quantities out of physical picture of the interactions. It means that di- 
mensionless quantities for the description of physical processes are used. The scaling 
function i]){z) depends in a self-similar manner on a single dimensionless variable z. It 
is expressed via the invariant cross section Ed 3 a /dp 3 as follows 

ns , T ^d 3 a 

Here, s is the center-of-mass collision energy squared, a in is the inelastic cross section, 
J is the corresponding Jacobian. The factor J is the known function of the kinematic 
variables, the momenta and masses of the colliding and produced particles. 
The function ip(z) is normalized as follows 



roc 

/ ^{z)dz = 1. 
J o 



(7) 



The relation allows us to interpret the function if)(z) as a probability density to produce 
a particle with the corresponding value of the variable z. 

We would like to emphasize that the existence of the function ip(z) depending on 
a single dimensionless variable z and revealing scaling properties is not evident in 
advance. Therefore the method proposed to construct ip(z) and z could be only proved 
a posteriori. 

2.3 Fractality 

Principle of fractality states that variables used in the description of the process 
diverge in terms of the resolution. This property is characteristic for the scaling variable 



The variable z has character of a fractal measure. For the given production process 
(1), its finite part zo is the ratio of the transverse energy released in the binary collision 
of constituents (2) and the average multiplicity density dN/dr] \ v= o. The divergent part 
Q^ 1 describes the resolution at which the collision of the constituents can be singled 
out of this process. The fl(xi,X2) represents relative number of all initial configurations 
containing the constituents which carry fractions X\ and X2 of the incoming momenta. 
The Si and S 2 are the anomalous fractal dimensions of the colliding objects (hadrons 
or nuclei). The momentum fractions x\ and x<i are determined in a way to minimize 
the resolution x-i) of the fractal measure z with respect to all possible sub- 

processes (2) subjected to the condition (3). The variable z was interpreted as a 
particle formation length. 

As we will show later the scaling function of high-pr particle production is described 
by the power law, i]){z) ~ z' 13 . Both quantities, ip and z, are scale dependent. There- 
fore we consider the high energy hadron-hadron interactions as interactions of fractals. 
In the asymptotic region the internal structure of particles, interactions of their con- 
stituents and mechanism of real particle formation manifest self-similarity over a wide 
scale range. 

2.4 Scale-relativity 

The properties of particle interactions in the space-time reflect symmetries of Na- 
ture. The principle of motion relativity has been used to formulate non-relativistic and 
relativistic theory. The special theory deals with only inertial coordinate systems while 
the expression of physical laws in the general theory of relativity should be written into 
any curvilinear coordinate system. Principle of general relativity states that "the laws 
of physics must be of such a nature that they apply to systems of references in any 
kind of motion" . Application of relativity principle can be extended to state of scale of 
reference system [13]. There are convincing evidence to consider that scale as well as 
other quantities characterizing a reference frame should be used to describe a particle 
state in a high-p T range. In this range elementary probes such as direct photons, high- 
er hadrons and jets are not point-like objects. They have a complicated structure. 
The last is resulted from interactions of quarks, gluons and heavy bosons which are 
fundamental objects of theory. Therefore experimentally measurable quantities should 



z = ZqVL 



(8) 



where 




(9) 



depend on a ratio between scales of a studied object and a probe. In other words the 
variables used in the description of the process depends on the resolution. 

A generalization of the motion-relativity principle to the scale-relativity principle 
requires that "the laws of physics must be of such a nature that they apply to systems 
of references in any kind of motion and whatever its state of scale" . 

In the framework of data ^-presentation the mechanism of particle formation is 
described by the scaling function ip(z). Both quantities, ip and z, depend on the 
resolution Q" 1 while the anomalous fractal dimension 5 as found from our analysis of 
numerous experimental data to be resolution independent. We consider that the study 
of ^-scaling of high-p r particle production over a wide kinematic range of and y/s 
and determination of 5 could give new insight in theory of scale-relativity [12, 13]. 

3 Properties of data ^-presentation 

In the section we present and discuss some properties of data pt- and ^-presentation 
of high-pr particle production in pp, pp and pA collisions. We show that the scaling 
functions for different processes reveal the same properties. They are the energy inde- 
pendence of ip, the power behavior of ip at high- ,2 and A-dependence. 

3.1 Energy independence of ip(z) 

It is well known that numerous experimental data on high-p^ particle spectra mani- 
fest the strong dependence on collision energy y/s. The effect enhances as the transverse 
momentum of produced particle increases. 

3.1.1 7r + -mesons 

Figures 1(a) and 2(a) show the dependence of the inclusive cross sections of n + 
and K + meson production in pp collisions on the transverse momentum p T at incident 
proton momentum p L = 70,200,300,400 and 800 GeV/c and an angle 9 cm ~ 90°. 
The data were obtained at Protvino [26] and Batavia [27, 28]. Transverse momenta 
of produced particles shown in Figures 1(a) and 2(a) change from 1 to 10 GeV/c. We 
would note that data [27] and [28] corresponding to the momentum pi = 400 GeV/c 
are complementary and are in a good agreement each other. As seen from Figures 1 (a) 
and 2(a) particle spectra have power behavior at pi = (200 — 800) GeV/c. The effect 
of kinematic boundary is visible at the end of spectrum at pi = 70 GeV/c. 

The energy independence of data ^-presentation means that the scaling function 
■ip(z) has the same shape for different y/s over a wide pr range. 

As seen from Figures 1(b) and 2(b) ^-presentation of the same data sets demon- 
strates the energy independence of ip(z) over a wide collision energy and transverse 
momentum range. We would like to emphasize that the data [28] used in our new 
analysis confirm our earlier results [5, 29]. 

3.1.2 7r°-mesons 

The PHENIX Collaboration published the new data [30] on inclusive spectrum of 
7r°-mesons produced in pp collisions in central rapidity range at RHIC energy y/s = 200 
GeV. The transverse momenta of 7r°-mesons were measured up to 13 GeV/c. 



The datapp- and ^-presentations for 7r°-meson spectra obtained at ISR [31]-[36] and 
RHIC [30] are shown in Figures 3(a) and 3(b). One can see that p^-spectra of 7r°-meson 
production reveal the properties similar to that found for charged hadrons. The new 
data [30] on 7r°-meson inclusive cross sections obtained at RHIC as seen from Figure 
3(b) are in a good agreement with our earlier results [10]. Thus we can conclude that 
available experimental data on high-p T 7r°-meson production in pp collisions confirm 
the property of the energy independence of tp(z) in ^-presentation. 

3.1.3 Charged hadrons 

The STAR Collaboration published new data [37] on inclusive cross sections of 
charged hadrons produced in pp collisions at RHIC energy ^J~s = 200 GeV. The RHIC 
data and other ones obtained at U70 [26], ISR [38] and Tevatron [27, 28] are shown 
in Figure 4(a). Charged hadron spectra were measured over a wide kinematic range 
y/s = (11.5 — 200) GeV and p T = (0.5 — 9.5) GeV/c. The strong energy dependence 
and the power behavior of particle ^-spectrum are found to be clearly. The energy 
independence of data ^-presentation shown in Figure 4(b) is confirmed. It is of interest 
to verify the asymptotic behavior of tp at y/s = 200 GeV and reach value of z up to 30 
and more. 

3.1.4 Direct 7 

Direct photons are considered as the best probes of constituent interactions at 
high-pT- The calculations of direct photon cross sections in the framework of QCD 
are developed in next-to- leading order [39]. The basic mechanisms of direct photon 
production are considered to be Compton scattering (gq — > jq), and annihilation 
process (qq — > 7(7). These are direct mechanisms of photon production. In high- 
Pt range direct photons are also produced indirectly via bremsstrahlung of quarks 
(qq — > qq'-f, qg — > qgj). The contribution of indirect mechanisms of high-pr photon 
production can be significant. However there are significant theoretical uncertainties 
due to the choice of structure and fragmentation functions and renormalization, fac- 
torization and fragmentation scales. Therefore any reliable estimates of direct photon 
cross sections allowing to reduce some of theoretical uncertainties is of interest. 

The results of data analysis for direct photon production in pp collisions are pre- 
sented in Figure 5. Inclusive cross sections versus the transverse momentum at \fs = 
(24 — 1800) GeV over the range pr = (4 — 110) GeV/c are shown in Figure 5(a). Data 
used in the analysis are obtained by UA1 [40], UA2 [41], UA6 [42], CDF [43] and DO 
[44] Collaborations. Data ^-presentation (see Figure 5(b)) demonstrates the energy in- 
dependence of the scaling function if) of high-p^ direct photon production over a wide 
kinematic range. 

3.1.5 Jets 

First observation of jets in pp collisions at SpS was considered as compelling confir- 
mation of parton structure of hadrons. In high energy collisions of hadrons copiously 
jet production due to hard parton scattering is observed. A jet represents a group of 
moving collimated particles. In the framework of QCD jets are distinguished to can 
be quark and gluon. Quark and gluon jets are initiated by fastest quark and gluon, 
respectively. It should be noted that the mechanism of jet formation is insufficient 



explored and not clearly understood till now. Therefore we hope that scaling features 
of jet production could be useful to obtain additional constraints for models of jet 
formation. 

In Figures 6(a) we show the invariant cross sections of inclusive jet production in 
pp collisions at y/s = 630 and 1800 GeV. These experimental data are obtained by 
CDF [45] and DO [46] Collaborations. A clear energy dependence of the cross section 
is observed to be. Difference between cross sections at y/s = 630 and 1800 GeV 
increases with transverse energy of jet. Data ^-presentation shown in Figure 6(b) 
demonstrates independence on collision energy y/s. The anomalous fractal dimensions 
for jet production in pp and pp collisions was found to be constant and equal to 1 [9]. 

3.2 A-dependence of ip(z) 

A comparison of particle yields in hadron-hadron and hadron-nucleus collisions is 
a basic method to study the nuclear matter influence on particle production. The 
elementary process is considered as probe of more complex system like nucleus. The 
difference between cross sections of particle production on free and bound nucleons 
was considered as an indication of unusual physics phenomena like EMC-effect, J/ip- 
suppression and Cronin effect. 

A change of the shape of px spectra is considered to be evidence that the mecha- 
nism of particle formation in nuclear matter is modified. Therefore it is convenient to 
compare scaling functions corresponding to different pA processes over a wide range of 
y/s and pr- 

The search for scaling features of particle formation in pA as well as in pp collisions 
and the study of their dependence on the atomic weight A is of interest for development 
of theory. 

A-dependence of ^-scaling of hadron production in pA collisions was studied in 
[7, 11]. It was established z-scaling for different nuclei (A = D — Pb) and types of 
produced particles (TT ±fi ) K ± ) p). To compare the scaling functions for different nuclei 
the symmetry transformation 

z — > ocaz, ip — > a^V (10) 

has been used. The parameter a of the scale transformation (10) depends on the atomic 
weight A. It was parameterized by the formula ct(A) = 0.9A 015 . 

Figures 7(a) demonstrates the spectra of 7r + -mesons produced in proton-nucleus 
collisions at y/s = 11.5 and 27.4 GeV and 9^J/ ~ 90°. Our new data analysis includes 
experimental data obtained at Protvino [26] and Batavia [27, 28]. The data [28] extend 
transverse momentum range up to px = 8.5 GeV/c. A good compatibility of [27] and 
[28] data sets in the overlapping region was observed. Solid and dashed lines are 
obtained by fitting of the data for W, Pb and D, respectively. They demonstrate the 
strong dependence of pr-spectra on collision energy y/s. 

Figure 7(b) shows the ^-presentations of the same data. The obtained results is 
the new confirmation of ^-scaling of high-pr hadron production in pA collisions. The 
universality of the scaling function ip for different nuclei means that mechanism of 
high-p T particle formation in nuclear matter reveals property of self-similarity. 

3.3 Power law 

One of the general properties of data ^-presentation is the power law of scaling 



function 

i;{z)~z- p . (11) 

Such behavior of ip as seen from Figures l(b)-7(b) is observed for different particles 
(hadrons, direct photons and jets) produced at z > 4. The data sets demonstrate 
a linear z-dependence of ip(z) on the log-log scale at high z. The quantity (3 is the 
slope parameter. The value of the slope parameter (5 is found to be constant with high 
accuracy and it is independent of energy y/s over a wide high transverse momentum 
range. Some indications (for tt° mesons, charged hadrons, direct photons, jets) are 
obtained that the value of slope parameter for pp is larger than the corresponding 
value for pp collisions, /3 PP > (3 PP . 

The existence of the power law, means, from our point of view, that the mechanism 
of particle formation reveals fractal behavior. 

4 Multiplicity charged particle density 

The important ingredient of ^-scaling concept is the multiplicity density of charged 
particles, dN/drj(s, rj). The scaling function ip and the scaling variable z is proportional 
to [dN/ drf\~ l . In the first case the quantity is included in the expression (5) to normalize 
the function ip and to give the physical meaning for it as a probability density to produce 
a particle with formation length z. In the second one (7) the multiplicity density is 
taken at i] = 0. Therefore z is proportional to energy of elementary subprocess per one 
particle produced in the initial hadron collision. 

The energy dependence of the multiplicity charged particle density for inelastic and 
non-single diffractive pp and pp collisions is shown in Figure 8(a) [47]. The collision 
energy sfs changes from 14 to 1800 GeV. New data for dN/df](s,i]) as well as for 
inclusive cross section Ed 3 a/dp 3 for pp collisions at RHIC energies are of interest for 
verification of z-scaling. 

5 7/V° ratio for pp and pp 

The properties of the scaling function for direct 7 and 7r°-meson can be used to 
estimate the dependence of the 7/71" ratio of inclusive cross sections on transverse 
momentum at LHC energies. 

The asymptotic behavior of ip(z) was found to be described by the power law for 
7r°-meson and direct photon production in pp and pp collisions. The slope parameters 
are satisfied to the relations fP m > (3^ p , f3 pp > /3p p , (3^ p > fP pp and j3^ p > fP pp . 

As seen from Figure 3(b) the cross section data [30] of 7r°-meson production in pp 
collisions obtained by the PHENIX collaboration at RHIC are in a good agreement 
with the asymptotic behavior of ip(z). 

Figure 8(b) shows the 7/71-° ratio of inclusive cross sections as a function of the 
transverse momentum pt at \fs = 5.5 and 14. TeV. The ratio was found to be different 
for pp and pp collisions. It increases with p?- The ratio has the cross-over point at 
p T ~ (60 — 70) GeV/c andj>r — (110 — 130) GeV/c for pp andpp collisions, respectively. 

6 z — Pt plot 

The z — Pt plot is the dependence of the variable z on the transverse momentum 



Pt of produced particle for a given process. The plot allows us to determine the high 
transverse momentum range interesting for verification of ^-scaling and experimentally 
inaccessible up to now. 

Figure 9(a) shows the z — pr plot for the pp — > n + X process at y/s = (24 — 14000) 
GeV. As seen from Figure 1(b) the scaling function tp(z) was measured up to z ~ 30. 
The function ip(z) demonstrates the power behavior at z > 4. Therefore the kinematic 
range z > 30 is of more preferable for experimental investigations of z-scaling violation. 
The boundary z = 30 corresponds to the different values of the transfers momentum 
Pt depending on collision energy y/s. 

Figure 9 (b) shows our predictions of the dependence of the inclusive cross section 
Ed 3 a /dp 3 on the transverse momentum p T for 7r + -mesons produced in pp collisions at 
the ISR, RHIC and LHC energies and an angle 9 cm = 90°. The verification of the 
predictions is of interest to determine the region of the scaling validity and search for 
new physics phenomena. 

7 "5- jump" 

Mechanism of particle formation at high transverse momenta in ^-presentation is 
described by the power law (11). Such behavior depends on the values of the anomalous 
fractal dimension of colliding particles, 8\ and 5 2 - The dimensions for hadrons, direct 
photons and jets produced in pp collisions were found to satisfy the relation 5h < <5 7 < 
5 j e t and to be independent of y/s and pr- The anomalous fractal dimension for nucleus 
8a is expressed via the dimension for nucleon 5n as follows 5a = A ■ Sn- 

Figure 10(a) shows the dependence of the anomalous fractal dimension 5 for 7r°- 
meson production in pp and pp [48] collisions on energy y/s. The value of 5h = 0.5 
used in our previous data analysis is confirmed by the new data [30] on inclusive cross 
section (see Figure 3(b)) obtained by the PHENIX Collaboration at RHIC. Figure 
10(b) gives evidence that values of the slope parameter j3 of the scaling function for pp 
and pp collisions differ each other at z > 6. 

The change of the fractal dimension 5 or "<5-jump" is considered as an indication on 
new mechanism of particle formation. It is assumed that the energy dependence of the 
quantity is especially sensitive in the high-pr range. Therefore the study of ^-scaling 
at higher y/s and px is of interest for search for new physics phenomena. 

8 Direct 7 and 77°-meson yields in pp and pPb 

collisions at RHIC and LHC 

The scaling properties of data ^-presentation for pp and pA collisions can be used to 
estimate particle yields in the kinematic region experimentally inaccessible at present 
time and to compare with other model predictions. 

Figures 11 and 12 shows our predictions of the dependence of the inclusive cross 
section Ed 3 a /dp 3 on the transverse momentum pt for direct photon (a) and ry°-mesons 
(b) produced in pp and pPb collisions at RHIC and LHC energies and an angle 9^/J/ = 
90°. The data on the cross sections [33, 42, 49] obtained at ISR energy y/s = (24 — 63) 
GeV are also shown for comparison. 



9 Conclusions 



The general concept of z-scaling for particle production in hadron-hadron and hadron- 
nucleus collisions with high transverse momenta was reviewed. The development of 
new method of data analysis based on data ^-presentation was presented. The scaling 
function ip(z) and scaling variable z were shown to be expressed via the experimental 
quantities, the invariant inclusive cross section Ed 3 a /dp 3 and the multiplicity charged 
particles density dN/dr)(s,r)). 

Physical interpretation of the scaling function i]){z) and variable z as a probability 
density to produce a particle with formation length z was argued. The quantity z was 
shown to reveal the property of the fractal measure and are the anomalous fractal 
dimensions of colliding particles. It was argued that ^-scaling reflects the fundamental 
symmetries such as locality, self-similarity, fractality and scale relativity. 

Results of new analysis of experimental data on the inclusive cross sections obtained 
at U70, ISR, SpS, Tevatron and RHIC were presented. The scaling properties of data 
^-presentation such as the energy independence, A-dependence and the power law 
were discussed. A complementary confirmation of z-scaling for 7r°-meson and charged 
hadron production in pp collisions at RHIC was obtained. 

New measurements of the multiplicity charged particle density and the inclusive 
cross sections of particle production in the experimentally inaccessible kinematic re- 
gion for the study of ^-scaling were suggested. The change of the anomalous fractal 
dimension ("5-jump") was suggested to consider as a new and complementary signa- 
ture of new physics phenomena of high-p T particle production in hadron-hadron and 
hadron-nucleus collisions at high energies. The z — pr plot was used to determine the 
regions that are of more preferable for experimental search for ^-scaling violation. The 
properties of data ^-presentation were used to predict high-pr particle spectra at RHIC 
and LHC energies. 

Acknowledgments 

The author would like to thank I.Zborovsky, Yu.Panebratsev, O.Rogachevski and 
D.Toivonen for collaboration and useful and stimulating discussions of the problem. 



References 

[1] J.D. Bjorken, Phys. Rev. 179 (1969) 1547; J.D. Bjorken, E.A. Paschanos, Phys. Rev. 
185 (1969) 1975. 

[2] R.P. Feynman, Photon-hadron Interaction (Benjamin, New York, 1972). 

[3] V.A. Matveev, R.M. Muradyan, A.N. Tavkhelidze, Part. Nuclei 2 (1971) 7; Lett. Nuovo 
Cimento 5 (1972) 907; Lett. Nuovo Cimento 7 (1973) 719. 

[4] S. Brodsky, G. Farrar, Phys. Rev. Lett. 31 (1973) 1153; Phys. Rev. Dll (1975) 1309. 

[5] I.Zborovsky, Yu.A.Panebratsev, M.V.Tokarev, G.P.Skoro, Phys. Rev. D54 (1996) 5548. 

[6] I. Zborovsky, M.V.Tokarev, Yu.A.Panebratsev, G.P.Skoro, Phys. Rev. C59 (1999) 2227. 

[7] M.V.Tokarev, I.Zborovsky, Yu.A.Panebratsev, G.P.Skoro, Int. J. Mod. Phys. A16 (2001) 
1281. 



[8] M.V. Tokarev, E.V. Potrebenikova, Comput. Phys. Commun. 117 (1999) 229. 

[9] M.V. Tokarev, T.G. Dedovich, Int. J. Mod. Phys. A15 (2000) 3495. 

[10] M.V. Tokarev, O.V. Rogachevski, T.G. Dedovich, J. Phys. G: Nucl. Part. Phys. 26 
(2000) 1671. 

[11] M.V. Tokarev, O.V. Rogachevski, T.G. Dedovich, JINR Preprint E2-2000-90, Dubna, 
2000. 

[12] I. Zborovsky, hep-ph/0101018; hep-ph/0206047. 

[13] L. Nottale, Fractal Space-Time and Microphysics. World Scientific Publishing Co.Pte. 
Ltd. 1993. 

[14] M. Tokarev, D.Toivonen, In: Proceedings of the XXXII International Symposium on 
Multiparticle Dynamics, Alushta, Cremea, Ukrine, 7-13 September, 2002, p. 158-162, 
eds. A.Sissakian, G.Kozlov, E.Kolganova, World Scientific Publishing. Co.Pre.Ltd. 

[15] M. Tokarev, G. Efimov, hep-ph/0209013. 

[16] P. Aurenche et al., hep-ph/0005114; S. Catani et al., hep-ph/0005025; S. Alekhin et al., 
hep-ph/0204316. 

[17] W.T. Giele, E.W.N. Glover, D.A. Kosower, Nucl. Phys. B403 (1993) 633; Phys. Rev. 
Lett. 73 (1994) 2019. 

[18] T. SjSstrand, Comput. Phys. Commun. 82 (1994) 74; hep-ph/0001032. 

[19] G. Marchesini et al., Comput. Phys. Commun. 67 (1992) 465; G. Corcella et al., JHEP 
0101 (2001) 010; hep-ph/0011363. 

[20] F.E. Paige, S.D. Protopopescu, BNL Report No. BNL38034 (1986); H. Baer et al., 
hep-ph/0001086. 

[21] J. Ranft, Z.Phys. C43 (1989) 439; A. Capella et al., Phys. Rept. 236 (1994) 225; J. 
Ranft, hep-ph/9911213; hep-ph/9911232. 

[22] B. Andersson, G. Gustafson, B. Nielsson-Almquist, Nucl. Phys. B281 (1987) 289; B. 
Andersson et al., Z.Phys. C57 (1993) 485. 

[23] X.-N. Wang, M. Gyulassy, Phys. Rev. D44 (1991) 3501; Comput. Phys. Commun. 83 
(1994) 307; nucl-th/9502021. 

[24] H.J. Drescher et al., Phys. Rept. 350 (2001) 93; F.M. Liu et al., hep-ph/0307204; K. 
Werner et al., hep-ph/0306151. 

[25] V.S. Stavinsky, Physics of Elementary Particles and Atomic Nuclei, 10 (1979) 949. 

[26] V.V. Abramov et al., Pizma v ZEFT 33 (1981) 304; Sov. J. Nucl. Phys. 31 (1980) 937; 
Sov. J. Nucl. Phys. 41 (1985) 700; Sov. J. Nucl. Phys. 41 357 (1985). 

[27] JW. Cronin et al., Phys. Rev. Dll (1975) 3105 ; D. Antreasyan et al., Phys. Rev. D19 
(1979) 764. 

[28] D. Jaffe et al., Phys. Rev. D40 (1989) 2777. 



[29] M.V. Tokarev, In: Proceedings of the 6 International Workshop " Relativistic Nuclear 
Physics: from Hundreds of MeV to TeV" (Varna, Bulgaria, September 10-16, 2001), 
El, 2-2001-290, Dubna, JINR, p.280-300; hep-ph/0111202. 

[30] H. Torii (PHENIX Collaboration), In: Proceedings of the 16th International Conference 
on Ultra- Relativistic Nucleus-Nucleus collisions "Quark Matter 2002", Nantes, France, 
18-24 July, 2002, p. 753-756, eds. H.Gutbrod, J.Aichelin, K.Werner, Elsevier Science 
B.V.; S.S. Adler et al., hep-ex/0304038. 

[31] A.L.S. Angelis et al., Phys. Lett. B79 (1978) 505. 

[32] C. Kourkoumelis et al., Phys. Lett. B83 (1979) 257. 

[33] C. Kourkoumelis et al., Phys. Lett. B84 (1979) 277. 

[34] C. Kourkoumelis et al., Z.Phys. 5 (1980) 95. 

[35] D. Lloyd Owen et al., Phys. Rev. Lett. 45 (1980) 89. 

[36] K. Eggert et al., Nucl. Phys. B98 (1975) 49. 

[37] J. Adams et al., hep-ex/0305015. 

[38] B. Alper et al., Nucl. Phys. B100 (1975) 237. 

[39] P. Aurenche et al., Eur. Phys. J. C9 (1999) 107. 

[40] C.Albajar et al., Phys.Lett. B209 (1988) 385. 

[41] J.A.Appel et al., Phys.Lett. B176 (1986) 239; J.Alitti et al., Phys.Lett. B263 (1991) 
544; R.Ansari et al., Z.Phys. C41 (1988) 395. 

[42] G. Ballocchi et al., Phys. Lett. B 317 (1993) 250; Yad. Fiz. 57 (1994) 1694; Phys. Lett. 
B 436 (1998) 222. 

[43] F.Abe et al., Phys. Rev. Lett. 68 (1992) 2734; Phys. Rev. D48 (1993) 2998; Phys. Rev. 
Lett. 73 (1994) 2662. 

[44] S.Abachi et al., FERMILAB-Pub-96/072-E; Phys. Rev. Lett. 77 (1996) 5011; S.Fahey, 
Ph.D. Thesis, Michigan State University, 1995; Y.-C.Liu, Ph.D. Thesis, Northwestern 
University, 1996. 

[45] F. Abe et al., Phys. Rev. Lett. 62 (1989) 613; 68 (1992) 1104; 70 (1993) 1376; 74 (1995) 
3439; 77 (1996) 438; 80 (1998) 3461. 

[46] D. Elvira, Ph.D. Thesis (Universodad de Buenos Aires, 1995). 

[47] W. Thome et al., Nucl. Phys. B129 (1977) 365. 
D.R. Ward, Report No. CERN-EP/87-178. 

F. Abe et al., Phys. Rev. D41(1990) 2330. 

M. Adamus et al., IHEP Preprint 88-121, Serpukhov, 1988. 

G. Alner et al., Z. Phys. C33 (1986) 2330. 

[48] M.Banner et al., Phys. Lett. B115 (1982) 59. 



[49] E. Anassontzis et al., Z. Phys. C13 (1982) 277. 



10 

1 tr 
-1 

10 

10" 2 
10" 3 
10" 4 
10" 5 
10" 6 

10" 7 

-a 

10 

10" 9 
10 1 
10 1 
10" 1 



-10 - 



p L (GeV/c) 1 

• 70 
■ 200 

« 300 

• 400 
i 400 

• 800 



P+P->7I*+X 

Ed 3 a/dp 3 
(mb/GevV 3 ) 



8 10 

p T (GeV/c) 




a) 



b) 



Figure 1. (a) The inclusive differential cross sections for 7r + -mesons produced in pp col- 
lisions at pi a b = 70, 200, 300, 400 and 800 GeV/c and 9 cm ~ 90° as functions of the transverse 
momentum pt- Experimental data are taken from [26, 27, 28]. (b) The corresponding scaling 
function ip(z). 
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Figure 2. (a) The inclusive differential cross sections for X + -mesons produced in pp col- 
lisions at pi a b = 70, 200, 300, 400 and 800 GeV/c and 9 cm = 90° as functions of the transverse 
momentum px- Experimental data are taken from [26, 27, 28]. (b) The corresponding scaling 
function ip{z). 
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Figure 3. (a) The dependence of the inclusive cross section of 7r°-meson production on 
the transverse momentum p± in pp collisions at y/s = 30, 53, 62 and 200 GeV and an angle 
cm of 90°. The experimental data are taken from [30]- [36]. (b) The corresponding scaling 
function ip(z). 
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Figure 4. (a) Inclusive cross sections of charged hadron production in pp collisions versus 
transverse momentum at U70, ISR, Tevatron and RHIC energies y/s = (11.5 — 200) GeV 
and an angle 9 cm of 90°. Experimental data are taken from [26, 27, 28, 37, 38]. (b) The 
corresponding scaling function ip(z). 
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Figure 5. (a) The dependence of the inclusive cross section of direct photon pro- 
duction on the transverse momentum pt in pp collisions at ^/s = (24 — 1800) GeV. 
Experimental data obtained by the UAl, UA2, UA6, CDF and DO Collaborations are 
taken from [40]- [44]. (b) The corresponding scaling functions. 
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Figure 6. (a) Inclusive cross sections of jet production in pp collisions versus 
transverse momentum at Tevatron energies ^/s = 630 and 1800 GeV and 6 cm ~ 90° 
obtained by CDF and DO Collaborations. Experimental data are taken from [45, 46]. 
(b) The corresponding scaling function ip(z). 
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Figure 7. (a) Inclusive differential cross section for 7r + -mesons produced in pA 
collisions at = 11.5 and 27.4 GeV/c and 6™ ~ 90° as a function of the transverse 
momentum p T . Solid and dashed lines are obtained by fitting of the data for W, Pb and 
.D, respectively. Experimental data are taken from [26, 27, 28]. (b) The corresponding 
scaling function ip(z). 
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Figure 8. (a) The multiplicity charged particle density dN/drj as a function of 
collision energy y/s at r\ = for and pp collisions. Experimental data are taken from 
[47]. (b) The 7/V ratio of inclusive cross sections versus the transverse momentum px 
in pp and pp collisions at y/s = 5.5 and 14. TeV. 
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Figure 9. (a) The z — p T plot and (b) the dependence of the inclusive cross 
section of 7r + -meson production in pp collisions on the transverse momentum px at 
v/i = (24 - 14000) GeV and an angle 6 cm of 90°. 
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Figure 10. (a) The dependence of the anomalous fractal dimension 5(s) on collision 
energy a/s. (b) The scaling function ip(z) of 7r°-meson production in pp and pp collisions 
on the transverse momentum p? & t energy ^/s = 30 — 200 and 540 GeV and an angle 
8 cm of 90°, respectively. The experimental data are taken from [30]- [36] and [48]. 
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Figure 11. The dependence of inclusive cross sections of direct photon (a) and 
?7 -meson (b) production on the transverse momentum in pp collisions at y/s = 
(24 - 14000) GeV. Experimental data are taken from [33, 42, 49]. Solid lines and 
points (*, A, *, x) are the calculated results. 
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Figure 12. (a) The dependence of inclusive cross sections of direct photon (a) 
and 77°-meson (b) production on the transverse momentum px in pPb collisions at 
y/s = (31 — 8800) GeV. Points (A, o, *,*, + ) are the calculated results. 



